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Abstract 

We consider the roughening of small angle grain boundaries consisting of arrays of dislocations 
and found two transitions, corresponding to fluctuations of the dislocations along and perpendicular 
to the boundaries. The latter contributes to a large scale fluctuation of the orientation of the crystal 
but the former does not. The transition temperatures of these transitions are very different, with the 
latter occuring at a much higher temperature. Order of magnitude estimates of these temperatures 
are consistent with recent experimental results from elasticity and X-ray measurements in solid 4 He. 
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Since the discovery of an increase in torsional oscillator frequency in solid 4 He at around 
200 mK|lj, there have been renewed interests in its low temperature physicalproperties. 

Much recent focus is on the role played by defects in this system, q, s, y, y, h q b For 

example, it is suggested that large angle grain boundaries can exhibit superfluid behavior % 
8[. X-ray measurements [2| found a change in the orientational fluctuation of the crystallites 
at around 1.75 K. Recently Day and Beamish[3] found a change in the shear modulus with 
the same temperature dependence as that for the decoupling in torsional oscillators. They 
ascribe this to a change in the mobility of dislocations. 

There usually is a finite density of dislocations and it is important to consider the 
long range elastic interaction between them. The simplest arrangement of a collection of 
dislocations comes from the small angle grain boundary (GB). These boundaries are pinned 
by the Peierls potential even in the absence of additional impurities. There has been much 
interest in the study of the pinning of an elastic two dimensional interface. Above the 
roughening transition temperature Tr the free energy to create a step becomes zero and the 
interface is depinned. This roughening transition temperature is a function of the strength 
of the pinning potential. Even as the strength of the pinning potential approaches zero, Tr 
remains finite. 

In this paper we examine if a GB can roughen and found that for a "electrically neutral" 
system, there are two roughening transitions, corresponding to the motion of the dislocations 
parallel and perpendicular to the boundary. The latter contributes to a large scale fluctuation 
of the orientation of the crystal as is observed in the X-ray experiments^- The other 
transition occurs at a much lower temperature, does not contribute to the large scale angular 
fluctuations but, because of the change in mobility of the dislocation, can cause a change 

n 

in the elastic coefficients, similar to that observed experiment ally [3]. Order of magnitude 
estimates of these temperatures are consistent with experimental results from elasticityQ] 
and X-ray measurements 3] in solid 4 He. We now describe our results in detail. 

A small angle GB consists of an array of dislocations with parallel Burger's vectors b. 
The trajectory of a dislocation can be represented by the positions of elements separated by 
lattice constants a z = a along it. We describe the configuration of the GB by the positions 
Cj of the elements of the dislocations. The elastic interaction between two elements of the 
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FIG. 1: Schematic diagram of the grain boundaries. 

dislocations is given by the formula 

V = «/(4tt) [b • b'/R + b Rb' • R/R 3 ] (1) 

where R is the distance between the elements. For simplicity we assume an elastically 
isotropic system. In terms of the Lame constants A and //, k — 4oq//(// + A)/(2/i + A). In 
addition, for each dislocation there is a core energy contribution V c = (E c /a )L that is 
proportional to its length L — J dl, where dl = (dc 2 x + dc 2 y + dc 2 ) 1 / 2 . Finally there is the 
pinning energy which we assume to be of a Peierl's form . U x cos(c X j/ao) + U y cos(c y j/ao). 

We take the dislocation line along z; the Burger's vectors along y. Thus the GB is 
in the xz plane. This is illustrated in fig. 1. The roughening transition depends on the 
form of the energy change for small devaitions Scj in the location of the dislocations. This 
deviations of the dialocations are along the x (parallel to the GB) and y (perpendicular to 
the GB) directions only. The fluctuation in the location of the dislocations Sr does not have 
a component along z. The energy change determines the "elastic" properties of the GB. 
We shall be interested in the Fourier transform 5c(q) = ^ - Scj exp(— iq ■ Cj)/y/~N. The GB 
is in the xz plane. Hence the wave vector q is in the xz plane only. The distance between 
the dislocations a x = a /9 is related to the misorientation angle 6 between the grains. We 
consider the simplest "electrically neutral" system with zero total Burger's vectors consisting 



of two grain boundaries of opposite Burger's vectors; one at y=0 and the other one at y=d. 
We focus on the lowest energy acoustic mode so that the displacements on the two boundaries 
are the same. 

A single grain boundary corresponds to two infinite crystals of different orientation 
joining together whereas a grain will be surrounded by boundaries of opposite orientations. 
We find that a single grain boundary does not roughen and thus focus on the pair, which 
also occur physically. 

The contribution to the energy change from the core energy is equal to 

5V C = J2^E c q 2 z [\6cM)\ 2 + IMq)| 2 ]. (2) 
q 

We next look at the contribution from V. From eq. (1) the change in V is given by 

8V = K/(4n)b 2 5c iq 5c h „ q {D s i3 + D' tJ ) (3) 

where the contribution from dislocations in the same GB (self) is 

DS a = ID 1 " cos (^ • ^Hv^-i/i? + v,vy/# 3 ] (4) 

R 

Similarly the contribution from the interaction energy between dislocations on different 
boundaries is given by 

D 'ij = - X} 1 ~ cos ^ ' ^)][V l V,l/i? / + ViVjU] (5) 

R 

U = (y + d) 2 /R' 3 , r' = (R 2 + d 2 ) 1 ^ 2 . D is very similar to the dynamical matrix for the 
two dimensional Wigner crystal which has been considered in detail with the Ewald sum 
technique by Bonsall and Maraduddin 9]. The two dimensional sums D s , D' can be evaluated 
in the same manner. 

We find in the long wavelength limit the x mode and the y mode are not coupled. For 
the "x" mode, the energy change is 

5V C + 5V* = {C x q 2 x + C z q 2 z )\5cM)\\ (6) 

C x = Kb 2 ^ 1 / 2 + 2)d/(a x a z ), C z = 0.5E C — 0.21K& 2 /(a z ir). Without the core energy contri- 
bution C z is less than zero and the lattice is unstable. This is consistent with the fact that 
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FIG. 2: The roughening temperature in units of Kelvin for the y mode (solid line) and the x mode 
(dashed line) as a function of the core energy normalized by koo- The experimental transition 
temperatures are indicated by the dotted line. 

a 2D rectangular Wigner lattice is unstable. Similarly for the "y" mode, the energy change 
is 

5V c + 5Vy=(C' xq l + C' zq l)\5c y ^ (7) 

where C' x = nb 2 ^ 1 ' 2 + 7/A)d/(a x a z ), C' z = 0.5E C + 3.6nb 2 /a z . The "elastic constants" C 
for the x mode are much smaller than that for the y mode. When the pinning energy are 
included, the energy we get is of the same form as that in the study of the roughening 
transition. 

For a single unpaied grain boundary, the elastic energy is of the first power in q. The cost 
of the long range fluctuation is higher and these fluctuations are supressed. As a result, an 
unpaired grain boundary does not roughen. We next estimate the roughening temperatures 
of our system. 

From the study of the roughening transition we find that, as the pinning Peierls potential 
approaches zero, the roughening transition is given by 

kT c = 2a z C/a x (8) 

where C is the geometric mean of the elastic constants. For example, for the x mode 
C = (C x C y ) l l 2 . The crystal symmetry of 4 He is HCP. In the present calculation, we have 
approximated it by an elastically isotropic system. Our goal is not to make an accurate 
prediction of the temperature but to find out if the physics discussed in this paper is of 
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relevance to the experimental system. We have used the following estimates of the Lame 



constants 



ll[ // = 0.72K/A 3 . A = 1.7K/A 3 . From these, we find nb 2 = 3.7 x 10 2 KA. 



We have used|l2j 9 = 0.2 degree. The linear spacing between the dislocations is thus 
°>x — a o/@ — 1029A. From an estimate of the dislocation density of 6 x 10 9 / cm 2 = l/(a x d), 
we estimate a mean spacing between the boundaries as d = 10 4 A/6. The core energy is not 



known. In two dimensional melting, it is suggested [13j that when E c > E c q = 0.056«:ao, the 
transition becomes first order. This provides a sense of scale for E c . In fig. 2 we show the 
roughening temperature for the x and the y modes (dashed and solid lines respectively) as 
a function of E c . When E c is around E c q, our estimate of Tr for the x mode is of the same 
order of magnitude as the experimental transition temperatures of Day and coworkers^. 
Tr for the y mode is of the same order as that of Burns and coworkers [2J and not a strong 
function of E c . The experimental results are indicated by the dotted lines. We next examine 
the physical implications of the transitions. 

We first show that only the upper y transition is connected with the large scale orientation 
fluctutaion of the transition. Our calculations is connected with the fluctuation of the 
position of the dislocations. We first relate these to the fluctuation of the atomic positions. 
The displacement u of the crystal at r caused by a dislocation located at position c can be 
written in Fourier transform as u(r — c) = J <ikexp[ik • (r — c)]u& where 

u(Jfe) = b[e y /k 2 - 0.5(1 - 2a)/(l - a)kk y /k\] 

a is the Poisson ratio. The change in the displacement as the dislocations on the two 
boundaries are moved from c to c + 5c contains contributions from both grain boundaries 
and is equal to 8u = 5ui(r) — 5u 2 (r); <5ui(r) = X]j[ u ( r — c j ~ ^ c j) ~ u ( r — c i)]> ^ u 2( r ) = 
^j[u(r — Cj — e y d — 5c j) — u(r — Cj — e y d)]. Cj is in the xy plane. 

Let us look at the angular fluctuation. We have 59 = (—d x 5u y + d y 5u x )/2. For the 
"parallel" wave vector k p in the xz plan, we are interested in J dy < (59)k p (59)-k p >= 
J dk y < \59k pt f Cy \ 2 > . It is straightforward to show that 5u(k) = u(k)[g(k) — g'{k)} where 
g — g' Pd k y d *Y1 G k • 5c^ p . We finally obtain, with a = x, y, 

f dy< (59 kp ) 2 >= F a {k) < \5c a (-k p - G)5c a (k p + G)\ > . (9) 

J a,G 

F a (k) = J dky[—k x u y (k) + k y u x (k)] 2 (k y d) 2 k 2 , G is a reciprocal lattice vector. Thus the angu- 
lar fluctuation of a grain can be related to fluctuation of the position of the GB. Subtituting 



in the expression for u(k), we get 

F a = b 2 k 2 x J dk y {k y d) 2 k 2 Jk\ 

Thus, for k z = 0, 

J dy < (56 kp ) 2 >^ kl < \5c x (-k x )5c x (k x )\ > (10) 

for the x mode 

J dy < (56 kp ) 2 >^ k 2 x n < \Sc y (-k x )5c y (k x )\ > /a (11) 

for the y mode. 

From previous calculation of the roughening transition [loj], in a purely relaxational model, 

<\5c(-k x )5c(k x )\>^ l/[A(k 2 x + C 2 )]- 

As the transition is approached, the "elastic" coefficient A and the relaxational rate are 
expected to exhibit square root cusps; with £ ps exp[c/(Tjj — T) 1 / 2 ] for T < Tr. 

In the roughened phase with £ _1 = 0, < \Sc(— k x )5c{k x )\ > l/(Ak x ). As k x approaches 
zero, f dy < {89u) 2 > becomes finite above the higher temperature roughened phase but 
remains zero between the two roughening transitions. This is consistent with our interper- 
tation that the higher rougheningtransition for the y mode corresponds to the observation 
from recent X-ray measurements 2[] Above the x mode roughening transition, because of the 
change in mobility of the dislocation, a change in the elastic coefficients results; similar to 
that observed experimentally[3|. 

Experimentally the addition of 3 He in ppm to ppb concentration changes the temperature 
dependence of the elastic modulus anomaly. The core energy can be changed. As we see in 
fig. 2, Tr can be a sensitive function of E c . The impurities can also change the effective 
magnitude of the pinning potential as follows. We represent the pinning potential between 
3 He atoms at s a and the dislocation line elements by V 1 = J2 a ^( c j ~ s ")- This appears 
as a factor [exp(—V l /kT)] av in the partition function where the square bracket means an 
impurity averaging over the positions s a . If we approximate this factor by an cumulant 
expansion, we get to lowest order a term of the form exp([V l ] av /kT). Because the 3 He 
appears as a substitional impurity, [V l ] av has the same periodicity as the lattice and just 
modifies the effective strength of the Peierl's potential U x>y . The next order in the cumulant 
expansion provides for an effective short range attraction between the dislocation lines, which 



can modify the elastic constants C and hence T R . This short range interaction is probably 
weaker than the long range elastic interaction in eq. (1). 

The main difference between the current system and the pinning of flux line lattices (FLL) 
in superconductors and incommensurate charge density wave (CDW) systems is the presence 
of the commensurate pinning potential in the current system so that the "ground state" is 
simple. For the CDW and the FLL system, the focus is on the nature of the ground state 
as a result of the competition between the random pinning potential and the elastic energy 
cost to distort the lattice. In those systems, the interest is usually on three dimensional 
objects whereas the current system is two dimensional. 

In conclusion we found that grain boundaries can roughen; there are two transitions. 
Estimate of their physical properties and their transition tempertaures are consistent with 
recent experimental observations. We do not completely understand what is the connec- 
tion between grain boundary roughening and supersolid behavior. We can think of many 
scenarios but further work is necessary to clarify if any of them is valid. 

STC thanks N. Mulders for helpful discussions. 
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